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E^ X E'^ superstring theories are anamoly free and hence very attractive [1,2]. Witten 
has shown that in such theories there is a model independent U{1) Peccei-Quinn symme- 
try. When this symmetry breaks one gets an axion with a decay constant related to the 
compactification scale because the non-renormalizable interactions of the model indepen- 
dent axion arise as a result of compactification [3]. Furthermore, there are several model 
dependent axion degrees of freedom. 

In E^ X E'^ models compactified on a Calabi-Yau manifold the model independent axion 
degree of freedom corresponds to -B^j/, (/U, z/ = 1, .., 4), the zero mode of the antisymmetric 
tensor field -Bma^ {M,N = 1, .., 10) which is essential for anamoly cancellation [3]. The 
model dependent axion is a linear combination of the zero modes Bmn (^, n = 5, .., 10) [4]. 
Neither of the corresponding PQ symmetries survive compactification and the above axions 
have decay constants related to the compactification scale [5]. 

If Eg breaks down to a non-abelian group there are two non-abelian groups today: 
SU{3)c and the non-abelian subgroup of Eg. The model independent and the model 
dependent axion degrees of freedom rearrange themselves to give two physical axions- the 
QCD axion a and the Eg axion a', with decay constants close to the compactification 
scale [6]. The energy density of the axions is proportional to A^Fa (a,/3 > 0), where 
A is the scale at which the underlying group gets strong. The QCD axion with Fa > 
10 GeV and A ~ 200 MeV will not have decayed by today and suffers an energy density 
problem. In this paper we study the cosmological and astrophysical constraints on the 
energy density of the Eg axion and on its decay to photons to give limits on A for the 
Eg axion. We find that A should be greater than 1.2 x 10 GeV. A is also the scale 
of gaugino condensation in the hidden sector which breaks supersymmetry [7] . Since the 
masses of supersymmetric particles in the observed sector are approximately 10" A /Mp^, 
this corresponds to supersymmetric particle masses greater than 1 GeV. The scale of 
supersymmetry breaking is given by Mg ~ (A /Mpi) ' . Thus, Mg should be greater than 
10 GeV. Our result agrees with A ~ 5 x 10 GeV obtained by setting TOgaugino — ''^W [''']• 
We find it significant that the result obtained by us using cosmological and astrophysical 
constraints agrees with the independent requirement from particle physics. 

In a companion paper [8], we set A equal to 5 x 10 GeV and study the dilution of 
the baryon asymmetry due to axion decays. We find that for this value of A, most models 
of baryogenesis can not tolerate the large dilution (by a factor of 10 ) of the baryon 
asymmetry. 

World sheet instanton effects may add terms to the potential of the model dependent 
axion [9,10] for most Calabi-Yau spaces (though not all). In such a case, the QCD axion 
will not get its mass purely from an FF term but can have large contributions from other 



non-derivative couplings. The Eg axion will be almost equivalent to ai below [11]. The 
mass of the Eg axion will probably not be greatly affected. Hence, we ignore these effects 
below. 

II 

In an Eg x Eg superstring theory compactified on a Calabi-Yau manifold and in which 
the Eg breaks down to a non-abelian group, the axion lagrangian in 4 dimensions is given 

by [6] 

where ai and a2 are the model independent and model dependent axion degrees of freedom 
and Fi and F2 are the corresponding decay constants. Both Fi and F2 are related to the 
compactification scale and are about 10 GeV [11]. F can be taken to be either the 
electromagnetic or weak or strong interactions field strength tensor. F' corresponds to 
hidden sector gauge field strength tensors. The lagrangian can be rewritten as 

L = lid^af + ^{d^a'f + ^^^aig'^FF) 

?2 _ 77l2 



+ (^^a') ig'^F'F' + g^\-^FF) (2a) 

V327r2F^ )V ^ F1+F2 ' 



where 



Fiai + F2a2 

(F2+F|)1/2^ 



(26) 



and 



/ _ F2ai - Fia2 
"" "(f2 + F|)1/2' ^'"^ 

F„ = i(F2 + F|)l/2, ^2d) 



F^ = FiF2/(f2 + f|)1/2. (2e) 



a is the QCD axion and a' is the Eg axion. a' is massless at high temperatures. But it 
acquires a potential and a mass m as the universe cools. The low temperature (T < A) 
potential for a' is calculated to be [6,12] 

V = ^Va^|1 -exp(ia7l5Fl)|l (3) 



Therefore the low temperature mass of the axion a' is 

2 T A" 

a' decays to two gluons, photons or Z's or to a VF pair through the couphng a' FF. 
The gluons, W and Z particles predominantly create jets of mesons and some baryons. 
At T ~ A, the non-abelian subgroup in the hidden sector gets strong leading to the 
hadronization of particles. The hidden sector hadrons will have masses of the order of A. 
Therefore we can ignore decays to hidden sector particles. A large fraction of the energy 
released as mesons is transferred to the photons, either through scattering or annihilations. 
We shall assume that only 20 per cent of the energy of the decaying axions is transferred 
to the photons i.e. / = 0.2. This is an extremely conservative estimate. 

The lifetime of the axion is given by 



r = 8200nH9t + 9^ + glr'^ = 5.2 x 10-l^s(,t + ,| + 4)'^/^ " f^ (5) 

^0 



TWA A^ 



where (71^2,3 are the couplings of t/(l)y, SU{2)i and SU{?,) respectively. We have set 
Ji.-Ji equal to 1. We have checked that induced decays to photons and gluons due to 
coherent axion field oscillations [13,14,15,16] are inhibited due to Hubble expansion. For r 
between 1 s and about 10 s, the axion mass ranges between 10 GeV and about 10^ GeV 
and the axion decays predominantly to mesons through gluons. For this range of lifetimes 
one can ignore gi and §2 and let (73 = 1 (or a^ ~ 0.1). For lifetimes greater than about 
10 s the axion mass is less than twice the pion mass and the axion can only decay to 
photons. If one rewrites eqn. (2a) in terms of Fp^F^Mi one gets the lifetime by replacing 
{gf + gl + glr^ m eqn. (5) by {Ae^K 

III 

When the Peccei-Quinn symmetry breaks at early times the axion field need not be at 
the minimum of its potential. At a temperature T^ the field will start to oscillate about the 
minimum of its potential with a period to" (T). (Above T^ the period of the oscillations 
is greater than the age of the universe H^ . At high temperatures the mass has a 
temperature dependence, just as for the QCD axion.) We study the zero momentum mode 
of the Eg axion field which can be treated as a condensate of zero momentum particles. 
(Higher momentum modes are redshifted.) The energy density of the axion field today is 
given by [13,14,17] 

Pao = mono = -mlAl (6) 
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no is the number density today and Aq is the amphtude of the oscillations today. A 
classical treatment is permitted as the number density is large [18]. The energy density 
can be rewritten in terms of parameters at T^ when the field starts oscillating as 

Pao = mon, {'^)$- lmom,Al (^) |. (7) 

Qs is the effective number of relativistic degrees of freedom used to calculate the entropy. 
The above expression assumes that the axion field loses energy only due to the expansion 
of the universe and the variation of the mass with temperature and not due to decays. 
(The astrophysical constaints on massive decaying particles is given in terms of limits on 
the ratio of the energy density of the massive particles today, if they had not decayed, to 
the photon number density today.) Tj is determined by the condition toj = 3H [13]. 



T 

-'-I 






1I/2 

(8) 



g^ is the effective number of relativistic degrees of freedom used to calculate the energy 
density. We shall take Qs^ = g^ c^nd gsj: = 3.9. By the time the axion field starts oscillating, 
the axion mass has already attained its low temperature value in (4). One can verify this 
by checking that m < 3H at T = A when the field attains its low temperature mass 
(Appendix B). In deriving (8) we have assumed that the universe is radiation dominated 
when the axion field starts to oscillate. One can verify this by checking that m > 3H at 
Teq when pa = Pr, i-G., the axion field starts oscillating before the axion energy density 
dominates the universe (Appendix C). Combining (7) and (8) gives 

Pao- 2 1/4 "^0 ^t 3/2- ^^> 

Ai is taken to be F^ ~ 10 GeV. In a non-inflationary universe, the spatial average of 
Aj_/F^ is 7r/v3. In an inflationary universe, setting A^ equal to F^ assumes that the 
universe lies in a typical post-inflationary region with Aj/F^ ~ 1. Linde has criticized 
this assumption. We refer the reader to Ref. [19] for a discussion of this issue. Since the 
number density of photons n^ = {1.2 /Tx)g^T^ 



Pa/n^ „ = 0.627r9/4 ^^ ^^ . (10) 



,A3/2^3/2 

This is the ratio of the axion energy density today in the absence of decays to the photon 
number density today. 



Our method will be to look at the astrophysical and cosmological constraints that 



apply to pa/n^ 



for different lifetimes r. For lifetimes between 1 s and 1.1 x 10 s we 

calculate the range of A from (5) and substitute that in (10) and compare with the limits 

on the energy density from observations. We show that axions with lifetimes in this range 

are ruled out by cosmological arguments for F'^ = 10 GeV. Larger lifetimes are ruled out 

by limits on the distortions of the microwave and gamma ray backgrounds and by closure 

arguments. They would also give too small masses for the gauginos and supersymmetric 

scalars. r less than 1 s corresponds to A greater than 1.2 x 10 GeV. 

IV 

The presence of the axion has different effects during different epochs in the evolution 
of the universe. The presence of axions during the neutron-proton freeze-out (~ 1 s) affects 
the ratio of neutrons to protons at freeze-out and thereby leads to a change in the amount 
of helium in the universe. Their presence during light element synthesis (~ 1 minute) 
affects the amount of deuterium and He. The decay of axions to photons or to other 
particles that transfer their energy to photons increases n^ and dilutes ry, the ratio of the 
baryon number density to the photon number density. This then alters the calculations 
of light element abundances. For 10 s < r < 10 s, the decay photons interact with 
the background thermal photons and create a cascade of photons with enough energy to 
photodissociate deuterium and helium. For lifetimes between 10 s and (re) combination 
the decay photons change the shape of the spectrum of the cosmic background radiation. 
Decays occuring after recombination change the diffuse extra-galactic gamma ray spec- 
trum. Decays of lifetimes greater than the age of the universe can lead to the axion energy 
density overdosing the universe. 

We now explore the astrophysical constraints on the energy density of the axions for 
different lifetimes of the axion. 

ls<r<10s: Increased expansion rate and dilution of rj 

The interaction rate for the weak reactions n + z/g *-^ p + e'^ and p -|- z/g *-^ n + e'^ 
goes as GpT while the expansion rate of the universe is H = [(87r/3)p/Mp,] '■/'', where 
p is the energy density of the universe [20] . When the weak reactions are faster than the 
expansion of the universe the ratio of the neutrons to protons is given by 



n 

— = exp 

P 



Am 

1^1 



(11) 



Am is the difference between the neutron and the proton mass. In the standard cosmol- 
ogy, at t ~ 1 s (T ~ IMeV), the reaction rates become less than the expansion rate of 
the universe and the neutrinos freeze out and the protons and neutrons are no longer in 
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chemical equilibrium. Thereafter the neutron-proton ratio does not change from its equi- 
librium value at the freeze-out temperature, except for free neutron decays. The precise 
temperature at which the freeze-out occurs is determined by the expansion rate of the 
universe which is influenced by the energy density of the axion. At slightly later times 
between 1 and 3 minutes (0.3 MeV > T > 0.1 MeV), the light elements deuterium, He, 
He and Li are produced. 

An increase in the expansion rate due to the axions causes the weak reactions above 
to freeze out at an earlier temperature. This leads to a a higher neutron to proton ratio 
and a higher percentage of He, if all the neutrons are converted into helium. However, 
the deuterium and He abundances are sensitive to the expansion rate at temperatures 
between 0.08 MeV and 0.04 MeV when the strong reactions that convert deuterium and 
^He into He freeze out. The amount of deuterium, ^He and He produced also depends 
on rjns, the baryon-to-photon ratio during nucleosynthesis [21]. For a given tjq today, the 
baryon-to-photon ratio before the axions decay is higher for higher axion energy density. 
The amount of He today increases for higher rjns while the amount of deuterium and He 
decreases because the rates of the reactions that convert deuterium and He into He are 
proportional to rjns- 

The constraint that He makes up less than 25 percent of the energy density of the 
baryons in the universe and that the ratio of the abundances of deuterium and He to the 
baryonic abundance is less than 1.0 x 10^ gives us limits on pa/n^ ■ For 1 s < r < 10 s, 

the limit from light element abundances implies that Pa/n^ is less than 1.5 x 10 /^ GeV 

to 8.4 X 10^^/^^ GeV [21]. But the axion energy density lies between 7.5 x 10^ GeV and 
1.6 X 10 GeV. Therefore this range of lifetimes corresponding to A between 1.2 x 10 GeV 
and 4.3 x 10^^ GeV and the axion mass between 1.2 x 10^ GeV and 5.5 x 10^ GeV is ruled 
out. 

10"^ s < r < 1.1 X 10 s: Photodissociation 

For decays that occur in the time interval 10 s < r < 10 s, the decay photons can 
initiate electromagnetic cascades that can photodissociate deuterium and helium [22] . The 
high energy decay photons interact with the thermal photons producing electron-positron 
pairs. The electrons and positrons then lose their energy by inverse-Compton scattering the 
thermal photons to higher energies. These photons will repeat the above process till their 
energies fall below E^ ~ ?Tig/22T, the threshhold for pair production, {me is the electron 
mass.) Photons with energy greater than E^ and Q, the threshhold for photodissociation 
of an element, are more likely to pair produce than interact with the nuclei of the light 
elements [22] . The spectrum of the electromagnetic cascade initiated by a decay photon 
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of energy E{> E*) is [23] 

dN^ fi^^^7^''^ ^orO<E^<E,/2 




dE^ ~ ) ^EEiE-^ for E,/2 < E^ < E, 

for E^ < E^, 



(12) 



where the cascade spectrum has been normahzed as 

The photodissociation of an element starts when E^^ becomes larger than the photodisso- 
ciation threshhold for that element. 

The mass fraction of an element is defined as X^ = nfgZj^/njsf, where Zj^ is the atomic 
number of the element, n^ is the number density of the element and njy is the number 
density of the nucleons. The change in the mass fraction is given by [22,23] 

dXk^_dna r dNv rE*(t) dN^XkOk^^ 
dt dt Jq dE VJq^ dEj Uecrc ^ 

ffjQ, dE^ ^e(rc J 

where <Jj—,k is the partial cross-section for the dissociation of element j to k, a^ is the total 
cross-section for photodissociation of element k and a(j is the Compton scattering cross- 
section on the thermal electrons with density Ue- Anticipating the stringent conditions on 
the energy density of axions we assume above that Ue is not changed significantly by the 
decays and cascade process. -^ includes the change in the number density of axion due 
to decays only (and not due to the expansion of the universe) . 



dt T \ T 

where n^(t) is the value of Ua at some very early time tj much less than r, scaled to the 
time of dissociation t. e^^'^ ~ 1. -^ is the spectrum of decay photons. To simplify 
our calculation we only include photons coming directly from the decay of the axions. 
-^ is then 26{E — mo/2). The ratio of the decay rates into photons and into gluons 
is ctj^j^/ag ~ 10^ . Therefore we include a factor of f^ = 10^ below. Including all 
the energy transferred to photons via all the axion decay channels can only tighten our 
constraints. 



Integrating (14) gives 



X;' 



n 



f- 



Pa 



7% 



where 



and 



/3fc(r) 



1 

TJt 



dtexp 






Qk 



E dE^ aciE-y) ^ 



Pj^kir) = ll^ dtexp(-^)|^ 



E^t) I dN^ (T^^k 



^ , ^ / 7T^ — dE^ . 

E dE^ crciEj) 



J- 



(16a) 



(166) 



(16c) 



i above represents the time when the photo dissociation starts, i.e., when the cascade cutoff 
energy E^ equals the photo dissociation threshold Q. tjq is the baryon-to-photon ratio today. 
Above we have used 

niit)/ne{t)=ni{to)/ne{to) (17) 



and 

7 7 

Mto) - g^iVo = g^0^7o- 

The subscript refers to quantities today and A^ to nucleons. 
t* is calculated from 



(18) 



f 



1.6 X 10~^hgs ^Z^- 



Mp(^ GeV 



m. 



1/4 



F^(m2/22g)3/2 



(19) 



where Q is 2.2 MeV for deuterium, 6.5 MeV for ^He and 20MeV for "^He [24]. We have 
used the time-temperature relationship obtained in eqn.(ll) of ref. [8] to obtain (19). gg 
is 3.9. The photodissociation of light elements only goes on till the energy of the decay 
photons itself falls below the threshhold for pair production. Then no cascade appears and 
photodissociation effects can be ignored. The time at which this occurs is calculated from 



where 



mo 



2 



1.2 X 10"^ 



me 
22T 



Ff GeV 



1/3 



(20a) 



(206) 



(20b) is calculated from (5). Combining (20) with the t — T relationship 



T = 



s 1.6 X 10-25 Mp/^ GeV 



9. 



1/2 



PI 



2/3 



(21) 



gives t = 1.1 X 10' s. For decays that occur after this time, photodissociation of hght 
elements is neghgible. 

For 10 s < r < 1.1 X 10 s, the hmits on Pa/n-j, come from the constraints on the 
increase in deuterium and He due to photodissociation of He. (The depletion of deu- 
terium and He is eclipsed by the increase in deuterium and He due to photodissociation 
of He as He is 10 times more abundant.) 



xf{D + ^He) - X\D + 3He) ^ -Y\^Re)^^^^^r^ 

7 rjo n^ 



h, (22) 





where r = (|(T^4Hc^n3Hc,p3H + ^ VHe^npi?)/VHe^aZ; - 0-5- 

13 are computed numerically in Ref. [23] for a radiation dominated universe. The f in 
our matter dominated universe are much earlier than those calculated in Ref. [23]. Hence 
/^^Hc ^'^ be greater in our case for lifetimes between 10 s and 10 s. We take /94jjg to 
be the peak value obtained in Ref. [23] of 1 GeV" . We take y*( He) < .25 which implies 
r\ < 1.3 X 10^^ and consequently X^{D + ^He) > 3.4 x 10^^. To be consistent with 
galactic chemical evolution {D + 3He)/H < 10"^ or X^ {D + ^He) < 2.3 x 10"^. Thus 
the upper limit on pajn^f is 1.0 x 10^ /^ GeV. But the calculated value of pa/n^f is 
between 1.6 x 10 GeV and 5.1 x 10 GeV. Hence the corresponding values of A between 
4.3 X 10 GeV and 2.0 x 10 GeV are ruled out. The corresponding values of the axion 
mass lie between 5.5 x 10 GeV and 56 GeV. 

r > 1.1 X 10^ s: 

For lifetimes greater than 1.1 x 10 s, the decay photons would distort the cosmic 
microwave background or show up in the gamma ray background spectrum. Such lifetimes 
would also gives supersymmetric masses less than 100 MeV. Lifetimes greater than the 
age of the universe would overdose the universe. Furthermore, if one assumes that the 
maximum baryon asymmetry rj = nj^/uj produced in the early universe is 0(1) [25,26] 
and that the baryon asymmetry today is 3 x 10^ , then constraints on the dilution of the 
universe imply that A must be greater than 7 x 10 GeV [8]. 

Conclusion: 

Thus we see that astrophysical and cosmological arguments imply that A is greater 
than 1.2 X 10 GeV or that Mg is greater than 10 GeV. Our result agrees with the value 
of 5 X 10 GeV for A obtained independently by setting supersymmetric scalar masses 
equal to m^Y. 
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Appendix A 

In this Appendix, we stiow that the universe is radiation dominated at the temperature 
T ~ A, when the axion attains it low temperature mass. 

At T = A, the radiation energy density is given by 

.2 

(Al) 











Prad 


^ A 4 




while the 


axion 


energy 


density is 


given 


by 












Pa = 


= mgn 


= 


\mlA^(K) 
900^^ m2; 


A2(A) 

i?/2 

^a 



{A.2) 

The low temperature mass toq used above is given by (4). g^ is 106.75. At T = A the field 
has not started oscillating (see Appendix B) and we assume that A{A) is approximately 
F^. Also, A < Mpi. Thus, pa < Prad at T = A. 

Appendix B 

In this Appendix we show that the axion field starts oscillating at a temperature Tj 
below T = A at which the axion attains its low temperature mass. We shall show below 
that m < 3H at T = A, i.e., the axion field has not yet started oscillating at T = A. 

Since the universe is still radiation dominated at T = A (see Appendix A), 

3H = 0.97r3/2^y2 A^ ^^^^ 

Mpi 

However, 

Since A < F^ and (7^ = 106.75, m < 3H at T = A. Thus the axion field starts oscillating 
at a temperature T^ < A, when H has decreased sufficiently to satisfy m = 3H. 
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Appendix C 

In this Appendix we show that the universe is still radiation dominated at the tem- 
perature Ti when the axion field starts oscillating. We actually show that the axion field 
is already oscillating by the time the temperature has decreased to Teq, when the axion 
energy density is equal to the radiation energy density. 

At T = Teq, 

Mpi 

Since A(Tcq) < A{Ti) = F^ < 0.2Mpi, m > 3H at Teq. This implies that the universe 
started oscillating at a temperature greater than Teq, when the universe was still radiation 
dominated. 
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